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We discuss the problem of spin-orbit interaction in a 2D chaotic or diffusive quantum dot in 
the presence of exchange correlations. Spin-orbit scattering breaks spin rotation invariance, and in 
the crossover regime between different symmetries of the spin-orbit coupling, the problem has no 
closed solution. A conventional choice of a many-particle basis in a numerical diagonalization is the 
set of Slater determinants built from the single-particle eigenstates of the one-body Hamiltonian 
,.0 ■ (including the spin-orbit terms). We develop a different approach based on the use of a good-spin 

D ' many-particle basis that is composed of the eigenstates of the universal Hamiltonian in the absence 

of spin-orbit scattering. We introduce a complete labelling of this good-spin basis and use angular 
£^ ■ momentum algebra to calculate in closed form the matrix elements of the spin-orbit interaction 

' in this basis. Spin properties, such as the ground-state spin distribution and the spin excitation 

function, are easily calculated in this basis. 
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The statistical properties of the single-particle levels and wave functions in a quantum dot whose single-electron 
dynamics is chaotic can be described by random matrix theory (RMT) [1, 2]. More precisely, the RMT description 
holds for ~ gx single-particle levels around the Fermi energy, when gx ~ \/n is the Thouless conductance of a dot 
with n electrons. In a dot that is strongly coupled to external leads (i.e., an open dot), the RMT description of 
non-interacting electrons can be used to derive the mesoscopic fluctuations of the conductance. 

However, if the dot is weakly coupled to the leads (i.e., an almost-isolated dot), the number of electrons in the dot 
i-^J | is quantized and it is necessary to take into account electron-electron interactions. The simplest model in that case is 
£^ ■ the constant interaction (CI) model, in which a classical charging energy e 2 n 2 /2C (C is the capacitance of the dot) is 
added to the single-particle Hamiltonian. This charging energy term is responsible for the observed Coulomb blockade 
, peaks in the conductance as a gate voltage on the dot is varied. In a chaotic dot, the heights of these conductance 

peaks exhibit mesoscopic fluctuations. The conductance peak height distributions in the absence and presence of a 
time-reversal symmetry-breaking orbital magnetic field were derived using RMT [3] and were found to be in overall 
J> ■ agreement with the experiments [4, 5]. However, experimental studies of the peak spacing statistics [6-9] and the 
| temperature dependence of the peak height statistics [10] suggested that it is necessary to include interaction effects 

■ beyond the CI model. 

A systematic approach to correlations in chaotic dots is based on a screened Coulomb interaction. The randomness 
of the single-particle wave functions induces randomness in the interaction matrix elements when the latter are 
\Q ' calculated in the single-particle eigenstates [11]. It is possible to separate these matrix elements into an average part 
and a fluctuating part. The fluctuating part is of the order 1/gr [12] and can be ignored in the limit of large gr- 

■ The average part survives in this limit and leads to the so-called universal Hamiltonian [13, 14]. While there exists a 
strong-coupling phase whose fluctuation properties differ from those of the universal Hamiltonian [15, 16], in this work 
we assume the weak-coupling phase in which the universal Hamiltonian description is valid. The interaction part of 
the universal Hamiltonian contains, in addition to the charging energy term (which characterizes the CI model), an 
exchange interaction proportional to S 2 , where S is the total spin of the dot, and, in the absence of a time-reversal 

q | symmetry breaking orbital magnetic field, a Cooper channel term. The latter is repulsive in quantum dots and can 
Q . be ignored. 

In the absence of spin-orbit interaction, the random matrix ensembles describing the single-electron Hamiltonian 
• i— i , are the Gaussian orthogonal ensemble (GOE) and the Gaussian unitary ensemble (GUE) in the absence or presence 
' of a time-reversal symmetry breaking orbital magnetic field, respectively. The exchange interaction modifies the 
H | conductance [17, 18] and affects significantly the statistical fluctuations of both the conductance peak spacings and 
- - > peak heights. Using a realistic value of the exchange interaction constant ( J s = 0.3A), good agreement was found [18] 
between theory and the experiments of Ref. [8] for the temperature dependence of the peak spacing fluctuations, 
explaining previous discrepancies with the predictions of the CI plus RMT model. Discrepancies between the CI 
plus RMT model and the experiment of Ref. [10] for the peak height statistics at low and intermediate temperatures 
(kT < 0.6A) are also explained by exchange correlations [18, 19]. 

The presence of spin-orbit coupling in a 2D quantum dot was shown to lead to new RMT symmetry limits [20]. 
In particular, there are two spin-orbit terms that modify the vector potential of the orbital field: a term a± that 
is proportional to the spin component s z (perpendicular to the plane of the dot), and a spin-flip term a||. In the 
absence of a field parallel to the plane of the dot, these two terms lead to two new symmetry classes, both for 
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conserved and broken time- reversal symmetry (i.e., in the absence and presence of an orbital magnetic field). The 
predictions of the modified single-electron theory were confirmed in experimental studies [21] of large open dots, in 
which electron-electron interactions can be ignored (spin-orbit effects are enhanced in large dots). 

The universal Hamiltonian of an almost-isolated dot and in particular the exchange interaction term are modified 
in these new symmetry limits of the spin-orbit interaction [22] . The mesoscopic fluctuations of the conductance peak 
heights and peak spacings in these limits were studied and found to exhibit an interesting interplay between the 
spin-orbit and exchange interactions. 

However, in the crossover regimes between different symmetries, the Hamiltonian of the dot is no longer universal. 
Such crossovers describe many-body problems that cannot be solved in closed form. For a numerical diagonalization 
of the Hamiltonian, it is necessary to calculate its matrix elements in a many-particle basis. A conventional choice of a 
many-particle basis is obtained by finding the single-particle eigenstates of the one-body part of the Hamiltonian and 
constructing a complete set of Slater determinants that correspond to the occupied single-particle levels. The one- 
body Hamiltonian is diagonal in this basis, and the main part of the calculation of the Hamiltonian matrix elements 
is the computation of the two-body interaction matrix elements in this basis. Such a basis was used in Ref. [23], 
which studied the combined effects of spin-orbit and exchange interactions on the (/-factor statistics in metallic 
nanoparticlcs [24] . In such 3D nanoparticles there are no intermediate symmetries and the only relevant crossover (in 
the absence of an orbital magnetic field) is between the GOE and the Gaussian symplcctic ensemble (GSE). A basis 
of Slater determinants was also used in studying the ground-state magnetization properties of a disordered quantum 
dot in the presence of exchange correlations and fluctuations of the off-diagonal interaction matrix elements [25] . 

Here we consider another choice of the many-particle basis. We use a good total spin basis composed of the 
eigenstates of the universal Hamiltonian in the absence of spin-orbit interaction. There arc several advantages in 
using such a basis. First, the calculation of the matrix elements of H requires the evaluation of the matrix elements 
of the spin-orbit terms, which are one-body operators, in contrast to the conventional basis in which it is necessary to 
calculate the matrix elements of the interaction, a two-body operator. Second, the spin-orbit terms can be expressed 
as tensor operators under spin rotations, and we can take advantage of angular momentum algebra in the calculation 
of their matrix elements. In particular, matrix elements of both the a_L and ay spin-orbit terms can be calculated from 
the same reduced matrix elements of certain tensor operators. Third, the eigenstates of the universal Hamiltonian 
can be organized in ascending energy and constitute a suitable basis for truncated calculations. In particular, an 
attractive exchange interaction brings down in energy high-spin states that in absence of exchange lie higher in energy. 
Finally, the good-spin basis is very convenient for calculating spin properties of the dot such as the ground-state spin 
distribution and the spin excitation function. 

The current work introduces the good-spin basis formalism. This approach is useful in various problems involving 
exchange correlations. For example, it would be interesting to study signatures of quantum critical fluctuations in 
the spin excitation function that are predicted in Ref. [26] for spin-orbit crossovers in the presence of correlations. 
We emphasize, however, that our formalism is quite general and is not limited to spin-orbit scattering, ft can, for 
example, be applied to a quantum dot in which the electrons interact ferromagnetically between themselves and 
antifcrromagnctically with a spin-1/2 Kondo impurity [27]. In this problem the local spin of the dot plays the role of 
the one-body operator (instead of the spin-orbit terms). 

II. HAMILTONIAN 

In the absence of spin-orbit interaction, the Hamiltonian of the quantum dot in the limit of large Thouless conduc- 
tance gx is given by the so-called universal Hamiltonian 

Here aj^ is the creation operator of an electron in an orbital level [i and spin cr/2 (a = ±1), and are spin-degenerate 
single-particle levels. The last term on the r.h.s. of (1) is an exchange interaction with S = \ ^ a g , aJ la <r aa >a tJ , a i 
being the total spin operator of the quantum dot (<r are Pauli matrices) and J s is the exchange interaction constant. 
The single-particle levels are spin-degenerate and correspond to the eigenvalues of a single-particle Hamiltonian 
h = (—ih'V — eA) 2 /2m + V(r), where A = Bs(es x r)/2c is a vector potential describing an orbital magnetic field B3 
perpendicular to the plane of the dot and V(r) is a confining potential. In a chaotic dot, these single-particle levels 
follow RMT statistics, i.e., GOE (GUE) in the absence (presence) of an orbital magnetic field. 

In the presence of spin-orbit interaction, the single-particle Hamiltonian is modified. Performing a gauge 
transformation to a locally-rotated spin frame [20], the single-particle Hamiltonian is found to have the form 
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(— iKV — eA — ax — ay ) 2 /2m + V(r), where 
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are two effective vector potentials. The term describes an effective spin-orbit magnetic field that has opposite 
directions for spin up and spin down electrons, while ay describes a spin-flip term. Ai, A2 are the spin-orbit scattering 
lengths for electrons moving along the principal crystallographic directions e l5 e 2 in the plane of the dot. 
The many-particle Hamiltonian (1) is modified to include the spin-orbit Hamiltonian H so 
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In a chaotic dot, the spin-orbit Hamiltonian H so can be represented in terms of random matrices of the appropriate 
symmetry. We have 



i? so = ia_ 
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where r'l =T X - iT 2 and T ± , T 1 , T 2 are real antisymmetric random matrices, and a±,Q!|| are coupling parameters 
associated with the aj_ and aii spin-orbit terms, respectively. 

Here we work in the basis of single-particle eigenstates /x in the absence of spin-orbit interaction. These are the 
eigenstates of a real symmetric matrix (GOE) or a complex hermitian matrix (GUE) in the absence or presence of 
an orbital magnetic field, respectively. If the variance of an off-diagonal element is chosen as a 2 / 2, then the mean 
single-particle level spacing in the middle of the spectrum is A = na/\ // 2N (N is the dimension of the random matrix). 
The random Gaussian matrices r x , Ti and T 2 are also chosen to have variance a 2 /2 for their off-diagonal elements. 

The strength of the spin-orbit coupling is characterized by two dimensionlcss crossover parameters 



x±_ = n 



a ± rms(r^) f - 
— - — = aaViv , 



(5) 



and 
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In the actual physical application, these crossover parameters are determined by the two spin-orbit energy scales 
e± and e« 
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Here Li, 2 are the corresponding linear lengths of the dot, A ~ L X L 2 is the dot's area, and k, k' are constants of order 
unity. 



III. GOOD SPIN MANY-PARTICLE BASIS 



A numerical diagonalization of a many-particle Hamiltonian of the type (3) requires a choice of a many-particle 
basis. In the conventional approach, the many-particle basis is chosen to be the set of Slater determinants of the 
single-particle eigenstates of the full one-body Hamiltonian Y], lr e^a^a^a + H so - However, here we use a good-spin 
basis determined by the eigenstates of the universal Hamiltonian (1) (in the absence of spin-orbit interaction). To 
find such a basis, we note that the occupation number operator fix = hx+ + n\— of any single particle orbital A is a 
scalar under rotations in spin space and therefore commutes with the total spin operator, [n\, S] — 0. We can then 



4 



characterize the eigenstates of the universal Hamiltonian by the orbital occupations n = {n\} (n^ = 0, 1,2), total 
spin quantum number S and the spin projection S z = M. Additional quantum numbers 7 are needed to distinguish 
between states with the same orbital occupations n and total spin S. 

The sequence of occupations n = {n\} can be alternatively specified by {A.;}{A;}, where Ai (A,) denote the 
corresponding numbers of singly-occupied (doubly-occupied) levels. The many-particle basis is then denoted by 

\njSM) = |{AJ{AJ 7 SA/) • (9) 



A. Labelling of the many-body states 



A complete labelling of the good-spin states requires the specification of the quantum numbers 7. These quantum 
numbers can be described by a tree graph, as we explain in the following. 

Consider a set of orbital occupations n with q singly occupied levels Ai, A2, . . . , X q where X q > ■ • ■ A2 > Ai. They 
describe a manifold of 2 q Slater determinants 

4xai°U-" a Ul{M»S = 0> (10) 

with Oi = ±1 (i = 1, . . . , q). Since the doubly-occupied levels Ai have spin zero, the construction of good-spin basis 
within this manifold is equivalent to the problem of coupling q spin-1/2 particles to total spin S. The total spin S can 
then take the values q/2,q/2 — 1, ...,q/2 — [q/2]. A unique labelling of the states requires us to choose a convention 
for the coupling sequence of the spin-1/2 particles. We assume that the spins are coupled in the ascending order of 
their orbital label 

Ai g> A 2 <8> • • ■ ® X q = ((■ • ■ ((Ai g> A 2 ) ® A 3 ) <&■■■) <g> X q ) . (11) 




FIG. 1: An example of a tree graph representing the normal-coupling scheme Ai ® A2 <8> A3 (g> A4 (g> A5 for A5 > ■ ■ • A2 > Ai 
singly-occupied levels. The total spin S of each state is indicated on the vertical axis. The d q (S) states with the same total 
spin S axe distinguished by the different paths 7 of intermediate spins. The internal points are shown by solid circles, whereas 
the terminal points are shown by open circles. 



Fig. 1 illustrates this coupling scheme in a diagram showing the total spin S versus the number of singly-occupied 
levels q. The diagram has a tree structure. We call the first single-particle orbital state Ai the root of the tree. The 
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root bifurcates into two points with total spin S = and S = 1. The point S = 1 bifurcates into two more points 
with S = 1/2 and S = 3/2, while the point S = can only connect to an S = 1/2 state. This process continues with 
each point bifurcating into two (except for the S = points) until we reach the states with q singly-occupied levels. 
The points on the tree with q singly-occupied levels are called terminal points (open circles) while all other points 
between the root and the terminal points are called internal points (solid circles). We observe that in general there 
are multiple states with the same spin S (and orbital occupations n). The multiplicity is given by the combinatorial 
number 

d ^ = (sU)-{s + Ul) • < 12 > 

These states are degenerate eigenstates of the universal Hamiltonian and we distinguish them by additional quantum 
numbers 7. A convenient choice is the path used to generate the state in the coupling sequence, i.e., the spin sequence 
of the internal points, 

( 7 , S) = (Si = 1/2, S 2 , S 3 , ... , S,-i, S q = S). (13) 

This set of quantum numbers corresponds to a complete set of commuting observables (S 2 2 , S 2 23 , •• • i S 2 2 ( q -iy S 2 ), 
where S\2...j are intermediate spin operators. 1 



IV. THE TRANSFORMATION FROM A NORMAL TREE TO A CANONICAL TREE 

We call the tree introduced in the previous section a normal tree. The orbital levels Ai, A2, • ■ ■ , X q are ordered in 
ascending order and the normal-tree states are given by 



|{Ai}7SM) x = |(((Ai®A 2 ) S2 ®A3) g3 ®---®A g ) , (14) 



s 

where for simplicity we have omitted the labels Ai of the doubly-occupied levels. 

However, calculations of the matrix elements of a^a^^r (which appear in the spin-orbit interaction) are best carried 
out in the canonical basis, i.e., the fi-tvee (or i/-trcc). This is the coupling scheme where the last coupled spin is that 
of the orbital /1 (u) corresponding to the orbital label of the creation (annihilation) operator. This last orbital /j, can 
generally be out of order, as is the case for Xj = fj, and j < q. 

The normal- and canonical-tree bases are related by a unitary transformation, which, as we now show, can be 
expressed in terms of the Wigner 6j symbols. The transformation between the normal and canonical bases can be 
constructed by successive applications of two basic operations: commutation and association. The commutation of 
two spins ji , J2 coupled to spin j'12 is simply a phase 

(— 1)* 

Kh® h)h2 -> 172® 31)3121 ( 15 ) 

where <j) = j± + j'2 — j'12, as is easily verified by the corresponding symmetry of the Clebsch-Gordan coefficients. The 
operation 1Z of association among three spins J1J2J3 coupled to total spin J 

((ji ® h) jl2 ® is) j ^ (ji ® (h ® j 3 ) J23 ) 7 (16) 

is given in terms of the Wigner 6j symbol 2 

K jl2 , j23 = (-i^+^+^ + ./ v / (2ji2 + 1)(2j23 + 1) jit h g j (17) 

The matrix 1Z, defined for given spins ji, j 2 , J3 and J, and labelled by the intermediate spin values ji 2 and j'23, is a 
real orthogonal transformation. 3 



1 We use Sj to denote the spin of the coupled spin angular momentum Si2...j = Si + S2 + ■ ■ • + Sj. 

2 Relation (17) essentially constitutes the definition of the 6j symbol. 

3 The inverse transformation ji(j2J3) — > (3132)33 is simply obtained by inverting the sign of the phase in (17). 
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We consider the unitary transformation which takes us from the normal basis to a canonical basis where the last 
coupled spin is /x = Xj (j < q), the rest of the spins being coupled in their normal order. Such a transformation 
can be obtained by successive applications of the exchange operation R which exchanges the coupling order of two 
neighboring spins 

((As^ ® A;) ® A i+1 ) -3- ((As^ <g>A i+1 )„, <g>Ai) , (18) 

where A$ i _ 1 denotes a sequence of orbital labels coupled to spin $s_i. The exchange R can be obtained by the 
following successive application of commutation, association and commutation 4 

((M-^) Sj A i+1 ) gj+i (A i+1 (A^A^)^ *C ((X^As^s, a/)^ (19) 

— * A. i+ i) s , Xi 

Using (15) and (17), we find 

R St>8l = Rs^s'iSi-uSi+i) = (-1) S >- 2S ^- S '^(2S> + l)(25i + 1) || || , (20) 

where R(Si-i, Si+i) is a 2 x 2 matrix which depends parametrically on two spins Si— i, As for any exchange 

operator, R 2 = 1. 

The transformation from the normal basis to a given canonical basis can be achieved by a succession of exchange 
operators. It is important however to keep in mind that the many-particle states are defined in terms of the product 
of fcrmionic creation operators and not simply spins. The series of exchanges only reorders the labels of the operators, 
but we also want to move the corresponding creation operator itself to the top of the string of operators. This results 
in an overall sign which is determined by the parity number (even or odd) of exchanges. To take into account this 
fermionic sign, we simply redefine R in (20) 

R SuS , ee Rs^Si-uS^) = (-l)*- 2 *+i-^+y(2S? + 1)(25 4 + 1) || |"J |J . (21) 

V. CALCULATION OF MATRIX ELEMENTS 
A. The spin-orbit Hamiltonian in terms of irreducible tensor operators 

The use of irreducible tensor operators simplifies the calculation of the matrix elements through the application of 
the Wigner-Eckart theorem. Since both at and a\ a = {— l) i - 1+rT ^ 2 a\- (J transform as irreducible tensor operators of 
rank 1/2 under spin rotations, we can form the following irreducible tensor-operators from their bilinear products 

(A^)l ee (at ® ~a v f m = 2(1/2, a/2 ; 1/2, </ /2 | S, m)a\ a ~a va . , (22) 

crrr' 

where £ = 0, m = or £ = 1, m = 0, ±1. Eq. (22) reads 

0V)o = ( a l+ a v+ + a l- a »-) ( 23 ) 

and 

0V)o = — = (a\ + Ov+ - al_a v -) ; (A^)\ = -a\ + a v - ; (A fll/ ) 1 _ 1 = a f u _a^ + . (24) 



4 There are other sequences of operations that lead to the same exchange, but they can be shown to lead to the same matrix through the 
use of Wigner 6j identities. The sequence given in Eq. (19) is the "shortest path." 
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Using relations (23) and (24), we can rewrite the spin-orbit interaction terms (4) as 

H so = iV2a ± ^T^(A^)l-lia ll J2 T lu(A li u)l + h.c.) . (25) 

Both the a±_ and ay spin-orbit terms are components of a E = 1 tensor and thus can change the total spin of the dot 
by AS = 0, ±1. However, the aj_ spin-orbit term is the m = component of this tensor, thus conserving S 2 , while 
the spin- flip an term is composed of the m = ±1 components, thus changing M by ±1. 



B. Calculation of the reduced matrix elements 



According to the Wigner-Eckart theorem, the matrix elements of a tensor operator in a good angular momentum 
basis can be factorized into a reduced matrix element and a Wigner 3j symbol 

\Suou\ _/ t\S'-m' I S' E S 



( 1 'S'M'\{A^ l \ 1 SM) = (-lf- M \^ M , ~ ^jh'S'WA^WjS). (26) 

The reduced matrix element (7'S' || A^ v || 7S) is independent of the magnetic quantum numbers M, M' and m. Once 
this reduced matrix element of A^ v is calculated, we can easily obtained the matrix elements of both the a± and a\\ 
spin-orbit terms, using (26) with m = and m = ±1, respectively (for E = 1). 

The reduced matrix elements of the spin-orbit operators A^ v can be grouped into four different classes depending 
on the initial occupations of the orbitals fj, (0 or 1) and v (1 or 2): 

(i) n v = 1 — ► n' v = , n M = — > n' = 1 

(m) n„ = 1 — > n'„ = , n M = 1 — > = 2 

(ra) n„ = 2 — > n„ = 1 , n M = 1 — > = 2 

(ro) n y = 2 -h. < = 1 , TV = ^ = 1 . (27) 

Although the operators A^ conserve the total number of particles n, they do not necessarily conserve q on the tree 
graph. While cases (i) and (iii) conserve q, case (ii) is characterized by a q — > (q — 2) transition and case (iv) by a 
q — > (g + 2) transition. In the following, we need to consider separately these four classes of transitions. 

Below we outline the calculation of the reduced matrix element for case (i) in detail, while for the remaining three 
cases we only present the results. 



1. Transition from, a singly- occupied level to an empty level 

The first case corresponds to the situation in which a particle is moved from a singly-occupied level v to an empty 
level fi. The evaluation of the matrix elements proceeds in two steps. First we consider the evaluation of the reduced 
matrix element (7'S' || A$ v \ \ 7S) C in the canonical basis. Then we transform to the normal tree by application of 
appropriate orthogonal transformations introduced in Section IV. Note that these transformations conserve the total 
spin, i.e. S = S and §' = S'. 

Consider the reduced matrix element (7'S' 1 1 A^ v \ \ jS) c where the initial and final states (as well as the intermediate 
states) are in the canonical basis. 5 The action of A^ 1 in the canonical basis is illustrated in Fig. 2. Depending on 
7, the spin S" of the intermediate state (after a particle is removed from the level v) can be either S" = S + 1/2 or 
S" = S — 1/2. We note that S" is uniquely determined given the path 7. In Appendix A we show that the canonical 
reduced matrix clement is given by 

(7'S' || A% || 7S) C = (_i)S'+s"-=+i/a V(2E + i)(2S + 1)(2S' + 1) ^ |,j . (28) 
For E = 1, the corresponding 6j symbol leads to the expected selection rule S' = S, S ± 1. 



5 all canonical matrix elements will be denoted by the subscript c or else by the orbital label of the last coupled spin. 
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Ai A2 A3 A4 A5 A6 A7 



FIG. 2: An example of a q — 7 tree, showing the action of 1 in the canonical basis. The thick solid line represents the 
initial state 7 with S = 3/2. Under the action of A^ 1 , this state is connected to states with S' = S (dotted-dashed path) and 
S' = S + 1 (dashed path). In the case shown, the intermediate spin in the transition is S" — S + 1/2 (determined by the initial 
path 7). 



To evaluate the reduced matrix element in the normal basis (j'S' 1 1 A^ v \ \ -fS), we expand the initial and final normal 
states in terms of the respective canonical basis 

\jSM) = ^U^SM) V (29) 

7 

\iS'M') = J2 V YyW S ' M '^> ( 3 °) 

1' 

where subscripts fi and v refer to the orbital which is coupled last. Here the transformation matrices U v and V 11 are 
independent of the magnetic quantum numbers M and M'. Then, using Eqs. (28), (29) and (30), we obtain 

h'S'\\A^ \\ 7 S) = (-lf+z+V V(2E + 1)(25 + 1)(25' + 1) 

xE(-D S "^V{l/2 1/2 ( 31 ) 

7,7' 

As is expected, no magnetic quantum numbers appear in (31). 

To calculate the unitary transformation U" explicitly, assume the spin we move out is at position i in the normal- 
ordered tree, i.e., A; = v. The transformation U v can then be viewed as a push-out operation 

|A X , A 2 , ...,\i = v,..., \ q -^SM) A T |Ai, A 2 , . . . , Aj_i, A i+X , . .., A g , \ = v;^SM) v . (32) 

Defining 7 = (S\ = 1/2, S2, S3, . . . ,S q = 5) and 7 = (Si = 1/2, 52, S3, . . . , S q = S), can be written in terms of 
the exchange operator defined in Eq. (20) 

U"^ = R Si! g.(Si^i, Si + i) Rs i+1 ,s i+1 \Si> Si+2) Rs i+ 2,Si +2 \Si+ii S i+3 ) ■ ■ ■ R Sci _ l g q _ l (S q _2,S q = S) . (33) 

The matrices R in (33) are 2x2 (because of the spin selection rules), and they depend parametrically on the 
neighboring spins. We should note that a special situation arises when the spin to be pushed out is at the root of the 
tree, i.e. i = 1. In that case, a simple commutation is required to move the spin to the position i = 2 and one can 
formally replace Rg.^Si-i, S i+ i) = (-1) S2 . 
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Similarly. can be viewed as a push-in operation 

|Ai, A 2 , . . . , X q = ft; fS"M% V |Ai, A 2j . . . , A,- = fi, . . . , A,, ; j'S'M') , (34) 

and is calculated to be 

V YY = - R s;_ 1 ,s;_ 1 (5 , 9_2, S' q = S') R s ^ ^ s , q 2 (S' q _ 3 , S'^i) Rs> q ^s> q _ 3 (S' q -4, S' q „ 2 ) 

■■■R^s'XS'j-i^'j+J. (35) 

Eqs. (33) and (35) allow us to represent the unitary push- in and push-out operations (which are of dimension d q (S) x 
d q (S)) in terms of a string of scalar operations. 




(a) (b) (c) 

FIG. 3: An example of spin-tree operations involved in the calculation of the reduced matrix element ("/' S' \ \ A^ 1 \ \ ■yS) for 
case (i) and S = S — 1. We consider q = 7 singly-occupied levels (the multiplicity of the states at the internal nodes k is not 
shown explicitly). In panel (a) the thick solid line represents the initial path 7. The paths 7 which result from the push-out 
operation of state v (at position 2 in the normal-ordered tree) are obtained by taking any combination of dashed lines that 
leads to spin S in the terminal node, starting from node 1. Note that there are nine such paths. In panel (b), we show by the 
thick solid line one of the possible paths 7 from panel (a) and consider the action of the operator A^ 1 . The action of the 
annihilation operator in A^ 1 can only result in a backward motion on path 7 (labelled by 1), and leads to intermediate spin 
S" = S — i. The action of the creation operator in A^ 1 can lead to both S' = S — 1 and S' = S. We choose the path (labelled 
by 2) that leads to the desired final state S' = S — 1. Finally, in panel (c) we consider the push-in operation to move /1 back to 
its position 3 in the normal-ordered representation (denoted by the second dotted vertical line). The paths generated through 
this procedure can be obtained by any combination of the dashed lines that starts from the terminal node at S and connects 
to the thick solid line at position 2. There are five such paths. 

Fig. 3 illustrates the spin-tree operations which enter the calculation of (7' 5' || A^ v || 7S). The initial path 7 = 
(Si, S2, ■ ■ ■ , S q = S) in a normal tree that connects the root to the desired terminal point (e.g., the thick solid line 
in Fig. 3(a)). The possible paths 7 = (Si, S2, ■ ■ ■ ,S q = S) which result when the level v (at \ = v) is pushed out 
to the terminal point to form a canonical ^-tree are represented by any combination of the dashed lines in Fig. 3(a). 
Note the simple rule that results from the form of the exchange operator R Sk g k (Sk-i, Sk+i)- For k < i, Sk = Sk- 
For k > i, the (two) possible values of Sk are determined by examining the values of Sk-i and Sk+i- The Wigner 6j 
selection rules dictate that Sk = Sk~i + \ f° r Sk+i = §k-i + 1, = Sk-i — \ for Sk+i = §k-i ~ 1 an d §k = Sk-i ± \ 
for S^+i = Sk-i- Thus, given the initial path 7, the set of paths 7 generated by the push-out operation is easily 
determined. These rules dictate that 7 and 7 are identical for k < i, that they may differ by at most AS* = ±1 in 
each internal node for k > i and that their final spins are identical. 

Next we consider the operation of A^ v in the canonical basis, connecting a canonical u- tree with terminal spin S 
to a canonical /i-tree with terminal spin S' = S, S ± 1 (the particular values of S' depend on In Fig. 3(b) we 

show by the thick solid line one of the allowed paths 7 generated in Fig. 3(a). For this path, the allowed values for S' 
are S, S — 1. Since we are interested in the final state S' = S — 1, we have S' = S — 1, 7' = (Si, . . . , S q -i, S — 1), and 
the canonical matrix element of A^ v is given by Eq. (28) with S" = S q -i and S' = S — 1. Finally, in Fig. 3(c) we 
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show by dashed lines the resulting possible paths 7' after using the push-in operation to transform from the canonical 
tree 7' (thick solid line) to the normal tree (for which Xj = fj,). 

Thus, starting from an initial state with quantum numbers (7, S), the non-zero matrix elements of A^ v in the 
normal basis can be found by choosing all sets of allowed paths 7, 7' and 7' (as described above). For each set, we 
attach the matrix element R Sk g k (Sk—i, Sk+i) to each internal node k > i, the matrix element Rgi s' k (^'k-i^ ^fc+i) 
to each internal node k > j, and Eq. (26) to the terminal node. For a given final state (p/',S'), we then sum over 
all intermediate paths 7 and 7'. In practice, the summation is carried out only over 7 (or cquivalently 7'). This is 
because 7' and 7 may differ only in the value of the terminal node S q which does not enter the summation. 

2. Transition from a singly- occupied level to a singly- occupied level 

In this case a particle is moved from a singly-occupied level v to a singly-occupied level [i. Contrary to the previous 
case, the number of singly-occupied levels changes from q to q — 2 (since level [i becomes doubly occupied), and the 
final state is represented by a shorter tree. Since both levels v and /i exist in the spin tree representing the initial 
state (7, S), there is only one canonical tree defined by A 9 _i = /i and X q = v. This canonical tree is obtained from 
the initial normal tree by applying the push-out operation on both singly-occupied levels v and /1. Since we want 
level v to be the last coupled orbital, we first perform the push-out of fi and then the push-out of v. This process is 
described by 

(7, S) ^ (7, S), % (f, SV ^ (7", S'% "-h (7', S') , (36) 
where 7" = (S[, . . . , S' x ) and S" = S' q _i. The corresponding reduced matrix element is then given by 

(l'S'\\A%\\ 7 S) = ]T U^ r ( 7 'S'\\A%\\iS% , (37) 

7,7' 

where U is the push-out operator defined in Eq. (33). The canonical reduced matrix element in this case is found to 
be 

( 7 >S' \\A% || iS') c = ^ (2E + 1)(25' + 1)(25" + 1) |J 2 ^ |,| . (38) 

3. Transition from a doubly- occupied level to a single- occupied level 

Here a particle is moved from a doubly-occupied level v to a singly-occupied level [i. The level exists in the initial 
tree (7, S) while the level v exists in the final tree (7', S'). We implement the following process 

(7, S) % (7", S'% £ (7, S) ^ (f, ^ (V, 5') , (39) 

where 7" = (Si, . . . , S q = S, S q ± 1/2), S" = S ± 1/2 and 7' = (S( , . . . , S' q , S' q+1 ). The corresponding reduced matrix 
element in the normal basis is given by 

( 7 'S' || A% || 7 5) = E ^"7^7'^'||^||75)c, (40) 
7". 7, 7' 

and the corresponding canonical reduced matrix clement is 

(7'SM|(A-) S || 7 S) C = (-1) S ^^ !/' 2 |/} ■ (41) 

^. Transition from a doubly-occupied level to an empty level 

Here a particle is moved from a doubly-occupied level v to an empty level fi. This generates a canonical i//i-tree 
with A g+1 = and A 9+2 = v - The normal tree of the final state is then obtained by push-in operations on both levels 
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v and [i (the final state has two more singly-occupied levels than the initial state). We thus implement the following 
sequence 

(7, 8) * ( 7 ", S") v (f, S'U ^' (7, S) v ^ (7', 5') , (42) 

where 7" = (Si, . . . , S q = S, S q ± 1/2), S" = S ± 1/2 and 7' = (Si, . . . , S q , S", S" ±1/2). It follows that the allowed 
values of S' are S, S ± 1 . The reduced matrix element can then be expressed in the following form 

( 1 'S'\\A%\\ 1 S)= V^V^(j'S'\\A^\\ 7 S) c , (43) 

7" ,7,7' 

and the corresponding canonical reduced matrix element is given by 

tf§'\\A%, \\jS) c = (-l) S '+ s - s +y (2S + 1)(2S' + 1)(25" + 1) ^ |,j . (44) 

VI. OTHER TECHNICAL ASPECTS 
A. Organization and truncation of many-particle states in the good-spin basis 

In practical calculations it is often necessary to truncate the many-particle Hilbcrt space. For that purpose the 
many-particle good-spin states are arranged in ascending order of their total energy in the universal Hamiltonian 
(including the exchange interaction). Below we describe briefly the implementation of an enumeration scheme which 
takes maximal advantage of the spin and fermionic structure of the many-body Hilbert space. 

The universal Hamiltonian eigenstates are labelled by \wySM). We first enumerate these states according to their 
"distance" from the Fermi sea as measured by the electron-hole excitation metric. Thus, 

\h 2 ; h; e 2 ; e) , (45) 

where h is a vector of occupied hole states, h,2 is a vector of doubly-occupied hole states, and e, e 2 describe similarly 
the electron states. Next we implement the spin-labelling scheme. Since the doubly-occupied levels have zero spin, the 
spin structure is completely determined the singly-occupied levels. This structure is best represented by a spin-tree 
7 = (Si, S2, ■ ■ • , S). The resulting good-spin eigenstates of the universal Hamiltonian are then arranged in ascending 
order of their energy E n ^s and truncation can be done by choosing a suitable energy cutoff. This scheme allows for 
inclusion of low-energy large-5 1 states (in the case of ferromagnetic coupling) already at the level of the unperturbed 
basis states. This is particularly advantageous near the Stoner instability, i.e., for J s ~ A. 

Convergence in a truncated basis can be improved by using banded matrices to describe the spin-orbit matrices T 1 - 
and rll in Eq. (4). In particular, we have observed better convergence for Gaussian-banded matrices. Banded random 
matrices with a band width of order \/N (N is the dimension of the matrix) have similar statistics as in standard 
RMT. 



B. Spin-dependent observables 

The advantage of using a good-spin basis is transparent in the computation of spin-dependent observables. Physical 
observables can be classified as tensor operators under spin rotations and it is sufficient to calculate their reduced 
matrix elements, irrespective of the magnetic quantum numbers. We have already used this in the calculation of the 
matrix elements of the spin-orbit interaction. 

Spin-orbit interaction breaks spin symmetry, and in particular the ground state does not have good spin. However, 
it is still of interest to find the ground-state spin distribution. To find this distribution it is necessary to expand the 
ground-state in a good-spin basis. In our approach, such an expansion is automatically obtained once the many-body 
Hamiltonian is diagonalized in the good-spin basis. 

Another interesting spin-dependent observable is the transverse spin correlator 

(46) 



S+(u>) = Im 

7T 



(S+(u)S-(-u)) 
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This quantity can be easily computed by expanding the ground-state |0) and excited states \I) in the good-spin basis. 
At zero temperature, we have 



VS{S + 1) - M(M + l)(I\njSM + l)(n-ySM\0) 

nfSM 



i 

where we have used the known matrix elements of S+ in the good-spin basis 



S(u - Ei) , 



(47) 



(n'^S'M'lS+lnjSM) = ^S(S + 1) - M(M + 1) <5 M ',a/+i<5s',s<V,A 



(48) 



VII. CONCLUSION 



The problem of spin-orbit scattering in the presence of exchange correlations does not have a closed solution in 
the crossover between different symmetries of the spin-orbit scattering problem. We have developed an approach to 
solve the problem numerically that is based on an unconventional choice of the many-particle basis. Rather than 
using a basis of Slater determinants built from the single-particle wave functions of the single-particle Hamiltonian, 
we use a correlated good-spin basis. This basis is composed of the eigenstates of the universal Hamiltonian which 
includes an exchange interaction in the absence of spin-orbit scattering. For a given set of occupation numbers of the 
spin-degenerate orbital levels, we construct many-particle states of good total spin by coupling the singly-occupied 
levels as spin-1/2 particles. To distinguish between states with same spin, we use spin-tree diagrams. 

The spin-orbit interaction terms are one-body operators, and they can be decomposed into sum of spherical tensors 
under spin rotations. We have used angular momentum algebra to calculate the reduced matrix elements of these 
tensor operators in closed form. In a basis that is canonical with respect to the transition operator, the reduced 
matrix elements can be expressed in terms of Wigner 6j symbols. To find the reduced matrix elements in the normal 
basis, it is necessary to make a unitary transformation from a normal tree to a canonical tree. The calculation of this 
unitary transformation involves a series of spin exchanges, and each such exchange can also be expressed in terms of 
a corresponding 6j symbol. 

The use of the good-spin basis has several advantages: 

• It is only necessary to calculate reduced matrix elements of the spin-orbit terms. Such reduced matrix elements 
are independent of the magnetic quantum numbers, and the actual matrix elements are easily calculated using 
the Wigner-Eckart theorem. Furthermore, both the and ay terms of the spin-orbit interaction are calculated 
from the same reduced matrix elements of the operators A^ 1 . The a^ and ay terms are described, respectively, 
by the m = and m = ±1 components of A^ 1 . The dependence of the (non-reduced) matrix elements on the 
specific m component of A^ 1 enters only through a corresponding Clebsch-Gordan coefficient. 

• The exchange correlation energy is fully included in the unperturbed basis of the universal Hamiltonian. Such 
a basis is particularly useful when the exchange coupling is ferromagnetic and near the Stoner instability. If an 
energy-cutoff truncation is necessary (e.g., in large model spaces), our truncated good-spin basis will include 
high spin states that lie low in energy because of the attractive exchange interaction. 

• Spin-dependent observables such as the ground-state spin distribution and spin excitation function are easily 
calculated since the calculated eigenstates are already expressed as a superposition of good-spin states. 

• Our approach is not limited to the spin-orbit coupling problem but is applicable for any perturbation that is 
added to the universal Hamiltonian. Any one-body operator can be decomposed into spin tensor operators A^ v 
with S = 0, 1, and the matrix elements of such operators have already been calculated here in a closed form. We 
can also express any two-body perturbation as a superposition of tensor products of two such one-body tensor 
operators and A^, v ,. 

In this work, we have presented the good-spin basis formalism and its technical aspects. This formalism is useful in 
the study of problems involving exchange correlations in quantum dots. One such problem is the spin-orbit scattering 
problem discussed in this manuscript. Ground-state spin distribution, spin excitation function and conductance peak 
statistics in almost-isolated dots can be studied as a function of the exchange constant J s and the scaled spin-orbit 
crossover parameters x± and xn. In particular, it would be interesting to compare with signatures of quantum critical 
fluctuations [26] that were calculated in the crossover between various symmetry limits of the spin-orbit scattering 
problem. 
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Appendix: Derivation of the reduced matrix element of A^ v in the canonical basis 

Here we derive the reduced matrix element ("f'S 1 \\ A? v \\ */S) c in the canonical basis for case (i), describing the 
transition from a singly-occupied level to an empty level (see Sec. VB1). To this end, we will compute the matrix 
element ^(7' S' M'\(A /11J )^ l \'-fSM) u . Using the expression for the irreducible tensor operator (A^)^ in Eq. (22) and 
the resolution of the identity, £ r , s „ M „ If S"M") (7" S"M"\ = 1, we have 

^'S'M'\(al®~a v fj^SM) v = £ (-l)W*(l/2, a/2 ; 1/2, -a'/2 | E, m) 

era' M" 

^S'M'\al^''S''M''){^'S''M''\a v ^SM) u , (49) 

where 7" = (Si = 1/2, S2, . . • , S q -i) is a normal tree, and 5"' ~ is the intermediate spin after a particle is 

removed from orbital v. Here, = n v = for \*/"S"M"), while n M = 0, n v — 1 for \"fSM) u and = 1, n v = for 
1 7' S' ' M') M . Since a^ a transforms as an irreducible tensor operator of rank 1/2 under spin rotations, we have 

^'S'M'\alJj"S"M") = (-l)*' 1 (S", M" ; 1/2, a/2 \ S'M>) , (50) 

and 

{i'S''M"\a va ^SM) v = v ^SM\al al \i'S"M")* = (—l) g ~ 1 (S",M" ; l/2,a'/2\SM) . (51) 
The factors (— 1) 9_1 account for the fermionic nature of the spin operators. Substituting these relations in Eq. (49) 

^S'M'\{al®~a v f m \ 1 SM) v = £ (-l)( W V 2 (l/2, a/2 ; 1/2, -a'/2 | E,m) 

00' M" 

x (S", M" ; 1/2, a/2 \ S'M')(S", M" ; 1/2, a'/2 \ SM) . (52) 
This can be rewritten in terms of Wigncr 3j symbols as 

MS'M'\ (at ® ~a v f m \>ySM)u = E (-l) s+s '- s "-^ M ' (-lf' + ^- M V(2S + 1)(25 + 1)(25' + 1) 

aa' M" 

x (i 2 a _!</ (cy m" _m) m') • ( 53 ) 

Using a sum rule relating 3j symbols to a 6j symbol [28] and symmetry properties of the 3j symbol, Eq. (53) can be 
expressed as 

vtfS'M'\ (at ® 5„)^ |75M), ^ (2S + i)(2 5 + 1)(25' + 1) 





' 5' 


£ S\ 




,S" 


H 






to My 


{f 


1 

2 


S"j 



(54) 



Comparing with Wigner-Eckart theorem, Eq. (26), we obtain Eq. (28) for the reduced matrix element. 
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